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We study the cylinder-plate and the cylinder-cylinder Casimir interaction in the {D + 1)- 
dimensional Minkowski spacetime due to the vacuum fluctuations of massless scalar fields. Different 
combinations of Dirichlet (D) and Neumann (N) boundary conditions are imposed on the two in¬ 
teracting objects. For the cylinder-cylinder interaction, we consider the case where one cylinder is 
inside the other, and the case where the two cylinders are outside each other. By computing the 
transition matrices of the objects and the translation matrices that relate different coordinate sys¬ 
tems, the explicit formulas for the Casimir interaction energies are derived. From these formulas, we 
compute the large separation and small separation asymptotic behaviors of the Casimir interaction. 

For the cylinder-plate interaction with R L, where R is the radius of the cylinder and L is the 
distance from the center of the cylinder to the plate, the order of decay of the Casimir interaction 
only depends on the boundary conditions imposed on the cylinder. The orders are L“^^^/ln(L) 
and L~^~^/ In L respectively for Dirichlet and Neumann boundary conditions on the cylinder. For 
two cylinders with radii Ri and R2 lying parallelly outside each other, the orders of decay of the 
Casimir interaction energies when R1 + R2 L are I/~^^^/(lnL)^, L“^“^/lnL and L~^~^ respec¬ 
tively for DD, DN/ND and NN boundary conditions, where L is the distance between the centers 
of the cylinders. The more interesting and important characteristic of Casimir interaction appears 
at small separation. Using perturbation technique, we compute the small separation asymptotic 
expansions of the Casimir interaction energies up to the next-to-leading order terms. The leading 
terms coincide with the respective results obtained using proximity force approximation, which is 
of order where d is the distance between the two objects. The results on the next-to- 

leading order terms are more interesting and important. We And some universal behaviors. It is 
also noticed that for the case of Dirichlet-Dirichlet cylinder-plate interaction, the next-to-leading 
order term agrees with that obtained using derivative expansion. Hence, based on our results on 
other boundary conditions and on the cylinder-cylinder interaction, we postulate a formula for the 
derivative expansion to expand the Casimir interaction energy up to the next-to-leading order terms 
for DD, DN, ND and NN boundary conditions, for the interaction between two curved surfaces in 
[D -\- l)-dimensional Minkowski spacetime. It is found that the postulate agrees with our previous 
results on the sphere-sphere interactions except when D = 4. 
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I. INTRODUCTION 

In the pioneering work [I|, Casimir proposed the existence of a force between two parallel perfectly conducting 
plates due to the vacuum fluctuations of electromagnetic fields. This gives rise to the concept of vacuum energy, 
which was referred to as Casimir energy. In subsequent years, Casimir effect has been generalized to any quantum 
fields, and it is a purely quantum effect. The idea of this Casimir energy is quite simple. The ground state energy of 
a quantum harmonic oscillator is not zero, but it is equal to Iiw/2, where w is the frequency. A quantum field can be 
considered as the superposition of an infinite number of quantum harmonic oscillators, each with a different ground 
state energy. Casimir defined the Casimir energy as the sum of the ground state energies: 

Eo«=i:^. (1) 
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This sum is divergent and regularization is required. However, in the existence of two objects (boundaries), one can 
obtain a finite Casimir interaction energy after subtracting away the Casimir self energies of each of the objects. 

The Casimir self energy of an object is of its own interest and it has been under active investigation (see for a 
review). On the one hand, it is closely related to the one-loop effective action and on the other hand, it has been 
proposed to be a candidate for the dark energy Since the advent of string theory, studying physics in higher 

dimensional spacetime has become a norm rather than an exception. There have been quite a number of works that 
explored the Casimir energies of rectangular cavities, spheres and cylinders in higher dimensional spacetimes 

Casimir effect is more interesting when there exist two interacting objects. In the last century, theoretical compu¬ 
tations of Casimir interaction were limited to the configuration of two parallel plates. However, the advancement in 
nanotechnology and Casimir experiments at the end of the last century have called for theoretical understanding of 
the Casimir interaction between any non-flat objects. About ten years ago, a major breakthrough in Casimir research 
was brought by a few groups of researchers |17 h 4^ , which have shed new light on the research of Casimir interaction 
between two objects. Using worldline numerics, multiple-scattering method or mode summation method, exact for¬ 
mulas for the Casimir interactions of cylinder-plate, cylinder-cylinder, sphere-plate, sphere-sphere configurations have 
been computed. This has enabled the more precise analytical and numerical studies of the nature and the strength 
of Casimir force. Nonetheless, these works have been limited to the (3-|-l)-dimensional Minkowski spacetime. 

Last year, we have taken the first step to understand the Casimir interactions between non-flat objects in higher 
dimensional spacetime (dll - li^l . We have considered the sphere-plate and sphere-sphere interactions due to the vacuum 
fluctuations of massless scalar fields in {D + l)-dimensional Minkowski spacetime and studied the dependence of the 
Casimir interaction on the dimension of spacetime. Since cylindrical objects played an equally important role in 
physics as spherical objects, we explore the Casimir interaction between a cylinder and a plate, and between two 
cylinders in {D + l)-dimensional Minkowski s pac etime in this work. 

We generalize the formalism established in [40| to compute the Casimir interaction energy between a cylinder and 
a plate, between two parallel cylinders where one lies inside the other, and between two parallel cylinders exterior 
to each other. We consider massless scalar field with combinations of Dirichlet (D) and Neumann (N) boundary 
conditions. The generic formula for the Casimir interaction energy between two objects can be written in the form 

h 

Ecs = — InTr (1 - , (2) 

Jo 

and is thus known as the TGTG formula. Here T* is related to the scattering matrix of object i and can be computed 
by matching the boundary conditions on the object. The matrix G*^ is the translation matrix that relates the wave 
functions of object i to the wave functions of object j. The nontrivial problem is to compute these T® and G®-^ matrices. 

After deriving the TGTG formulas for the cylinder-plate and the cylinder-cylinder interactions, we derive the large 
separation and small separation asymptotic behaviors of the Gasimir interactions. The large separation asymptotic 
behavior is easy to compute since it only depends on a few entries in each of the T and G matrices. To compute the 
small separation asymptotic behavior beyond the leading term is a tedious task [s^, IdTI - f^ . As a confirmation of 
the correctness of the TGTG formula, the leading term of the small separation expansion of the Gasimir interaction 
energy is found to agree with that derive using proximity force approximation. One of the major contributions of 
the present work is the result of the next-to-leading order term of the small separation expansion. For the cylinder- 
plate configuration with DD boundary conditions, we find that our result agrees with that computed using derivative 
expansion in (s^ . Inspired by the work (53 |. we use our results on the cylinder-cylinder interaction to postulate a 
derivative expansion formula for the Gasimir interaction energy in (D + l)-dimensional spacetime, up to the second 
order term, for the interaction between any two objects with combinations of Dirichlet and Neumann boundary 
conditions. This ansatz is found to agree with the results we derive for the sphere-sphere interaction in except 
when D = i. 

This work will be interesting to those that wish to understand quantum field theory in higher dimensional spacetime. 


II. THE CASIMIR INTERACTION ENERGY 

In this work, we consider the vacuum fluctuations of a massless scalar field in {D + l)-dimensional Minkowski 
spacetime with metric 

ds^ = dt^ — dx\ — ... — dx\,. 

The equation of motion of the scalar field </j(x)e“®‘®®‘, x = (xi,..., xd) is 

/ 9 ^ 9 ^ \ 


( 3 ) 
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We will consider the following three problems: 

• The Casimir interaction between a cylinder and a plate. 

• The Casimir interaction between two parallel cylinders, one is inside the other. 

• The Casimir interaction between two parallel cylinders exterior to each other. 

The boundary conditions on the cylinder and the plate are either the Dirichlet boundary condition cj L , = 0 or 

.y X v' r I boundary 

the Neumann boundary condition = 0, where n is the unit vector normal to the boundary. 

We will take a cylinder to be 

2 I 2 _ 7d2 

Xi X2 = n , 

where R is the radius of the cylinder. Therefore, it will be convenient to work with the cylindrical coordinates where 


xi=rcos0, X2=rsin0, (4) 

so that the cylinder ccf + is given by r = i?. In the cylindrical coordinates, the equation of motion ([3]) reads 

as 


1 d 

dr'^ r dr 


1 92 Q2 q2 \ 



( 5 ) 


Solving this equation of motion, we find that the cylindrical waves can be parametrized by (n, fca,..., = (n, fcj_), 

where n is an integer, k± = (fca,..., kn) G and they can also be divided into regular and outgoing waves. The 

explicit formulas for these cylindrical waves are given by 


<fc,(x)=C:z:(Ar)e 


in6-\-ik^x^-\-.. ]yxD 


( 6 ) 


where * = reg or out for regular or outgoing waves, 


are normalization constants, and 




^OUt _ Z^g^+1 




A. The Casimir interaction energy between a cylinder and a plane 

For the Casimir interaction between a cylinder and a plate, we will take the cylinder to be 

Xi+ X 2 = R^, - - < Xi < —- for 3 < i < D, 

and the plate to be 

xi = L, —^ <Xi< for 2 <i < D. 

Here L > R and d = L — R is the distance between the cylinder and the plate. 

For the plane xi = L, we will parametrize the plane waves by the momenta perpendicular to the plane 
(^2,^3, ■■■,kD) = {k 2 ,kj_^). Solving the equation of motion ([3|) give the plane wave basis 

_ —isgn^kiXi+ik2X2+ik3X3 + ...+ikoXD ('7\ 

'^/C2,fc±V^l — 6 1 1C 


ki 



— 1-2 _ _ ^2 
t%2 ... ^r), 


where 
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and 


SgHres = 


SgHout = -1- 


Let H = H 3 ... Hjj. In the region between the cylinder and the plate, we can write the scalar field i^(x, t) in terms 
of the cylindrical coordinate system centered at the origin: 




n— — oo 


or in terms of the rectangular coordinate system centered at O' = Lei: 

dkjy 




( 8 ) 


(9) 


Here x' = x — Lei. 

Using the representation (|S]), we find that the boundary condition on the cylinder r = R gives rise to a relation of 
the form 


_ _ rpn^n,k± 

For Dirichlet (D) and Neumann (N) boundary conditions, T” is given by 


rpn,D,.^\ Ini'jR) 
Ti [lO =- 




( 10 ) 


( 11 ) 


K = 


KnijR)’ 

respectively. Here ^ is the imaginary frequency so that = uj, 

7 = \J+ kl + ... + kjj, 

Under the representation (jO)), the boundary condition on the plate x'l = 0 gives rise to a relation of the form 

^ _£fc2jfc2.fc_L_ (12) 

For Dirichlet (D) and Neumann (N) boundary conditions, is given by 

=1, = -1 (13) 


respectively. 

The two representations (|H]) and ([5]) are related by translation matrices V and W: 


n— — QO 

/ oo 77 


(14) 


It is easy to see that the matrices V and W are diagonal in kj^. In fact, canceling out fj-om both 

sides, we obtain exactly the same equation as in D = 3 dimension. Hence, quoting the result from dimension D — 3 
(see for example [1^), we have 


^n,k 2 
11^/C2 ,n 


( \/7^ + fcj + fc2 \ -^ 72 +fe|L 

^ ( V7^ + kl + k2 \ 
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Notice that (IT4l) implies that 




dk 2 

27r 




fe2,fcj 


^ Wk,,nV 


n,fe I 


n=—oo 


From (jlOF (1121) and (USD, we find that 

yn,k± _ _ rpn^n,k± 


dk 2 , 


= -T^H 2 I 

J — oo 


=T^H 2 


27r 




--T^H 2 r 

n' = —OO 


This is a relation of the form 


(15) 


where 


(I-M)B = 0, 
M = T 1 VT 2 W, 


(16) 


and B is the column matrix with components . 

The matrix B must be a nontrivial solution of ()16p . Hence, we obtain the dispersion relation 

det (I - M) = 0. 

Using standard contour integration technique, and the fact that all the matrices Ti, V, T 2 and W are diagonal in k±, 
we then find that the Casimir interaction energy is given by 

hcH dks dkj 

L ■./ 


ddcas — 


2 tt 


Trln(l 


where 


=T^H 2 


nrpk 2 


^^2 -T r rTik2 TT/ 

2^ Ui,fc2^2 ^k2,r. 


= T^Tr 


dk 2 


n-\-n' 


^ 7-00 2-^7^ + 
=T^f!^^K^+^,{ 2 ^L). 


\/ 7 ^ + fc| + fc2 \ -2^72+felL 


Recall that 

7 = + k\, + ... + fc|,. 

Since the dependence of Mn,n' on (fca,..., fc^i) is only through fc_L, we have 


ddcas —■ 


HcH 2 tt 2 r 

^ { 2 n)D- 2 r (i2^) X 

HcH 


dK dk± k^-^Ti In (l-M) 


2 D-2^§Y (-2^) 

HcH 


dk± k 


D-3 


D 

foo /.7 


d'y- 


d'y^ / d/cj_ 


fe 7 / 7 ^ - k']_ 

Trln(l - 


Trln(l-: 


D-3 

L 


2^ 2^°r(-2^)7o 7o 7/72 - 


2D- 


(17) 



















6 


B. The Casimir interaction energy of one cylinder inside the other 

For the Casimir interaction between two cylinders, one inside the other, we take the smaller cylinder to be 

2 I 2 7~)2 

Xi+ X 2 = rii, 

whose center is at the origin O, and the larger cylinder is taken to be 

{xi — LY + x\ = i?2, 


H, Hi 

—^ < Xi < — for 3 < i < D, 


< Xi < ^ for 3 <i < D, 


whose center is at O' = Lei. Notice that L < R 2 — R 1 and d = R 2 — Ri— L is the distance between the two cylinders. 

As in the cylinder-plate case, we can represent the scalar field (/^(x, t) in the region between the two cylinders in 
two different ways, one with respect to the cylindrical coordinate system centered at O; 




n= —00 


and one is in terms of the cylindrical coordinate system centered at O' = Lei: 


' f'i H ( "Aj- reg , /N out / /A -ruit 


n— — oo 


The boundary conditions on the cylinders give 

with 

InARl) 

H-n \lRl ) 

;pn,D/.>N Ar„(7i?2) 

T 2 (*4) = J , p X , 

InhR2) 

The two representations (fT51) and (ITOl) are related by 


^n,k^ = -Tn"d" 


Ti’A^O = 


AARi) 

KAlRiV 

KAR 2 ) 

IA 1 R 2 ) ■ 


AAkA^')= E 




AkAA= E 


n — — 00 


Compare to the D = 3 case (see e.g. I40|), we find that 

Vn,n' = Wn',n = (-I)”"”'(7A) • 

As in the cylinder-plate case, we then find that the Casimir interaction energy is given by 

hcH 


HCas — 




rm 


nOO 

/ d77^-2Trln(l-M) 

Jo 


(18) 


(19) 


( 20 ) 


( 21 ) 


( 22 ) 


(23) 


(24) 


where 


00 

Mn,n'=T^ E Vn,n"fA'Wn",u' 

n" — — oo 
00 

TE // 

T 'Y'n T 
^n" —n ^2 ^n"—n'’ 

n" — — oo 


(25) 
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C. The Casimir interaction energy of two parallel cylinder outside each other 

For the Casimir interaction between two parallel cylinders exterior to each other, we take one cylinder to be 

2 I 2 7t2 

Xi+ X2 = -ft^, 

whose center is at the origin O, and the second cylinder is taken to be 


H, Hi 

—< Xi < — for S <i < D, 


(xi — + xl = R? 


'2 I 


- - < Xi < —- for 3 < i < D, 

2 ~ 2 - - > 


whose center is at O' = Lei ■ In this case, L > R 1 +R 2 and d = L — Ri — R 2 is the distance between the two cylinders. 

In the region between the two cylinders, the scalar field (/3(x, t) can be represented by p8l) using the cylindrical 
coordinate system centered at O, or by (flUl) using the cylindrical coordinate system centered at O'. 

The boundary conditions on the cylinders give 


with 




'T'Ti n,k± 

ii a , 


(26) 

InilRi) 

~Kr^{-^R^y 


(27) 


In the present case, the two representations (ITS)) and m are related by 

00 

n= —00 

n' = —00 

In other words, 

00 

n' = —00 
00 

n— 

Compare to the D = 3 case (see e.g. d^), we find that 


n— — <X) 


(28) 


(29) 


O = (-1)” K^-n'hL). 


(30) 


From (l26)) and (|29ll . we find that 


j^n,k± _ _ rj^n^n,k± 

00 

= -rr ^ 


n — — CO 


L 2 C 


=T^ E 

n" — — co 

00 00 

=Tr E Kn"Tf E 


n — — CO 


n — — CO 


Then as in the cylinder-plate case, we find that the Casimir interaction energy is given by 

hcH 

TRZTT 


Ecas — 




pCO 

/ d77^-2Trln(l-M) 

Jo 


(31) 





where 


Mn,n'=Tr 

n" — — oo 
oo 

=T” ^ Kn-n"Tf K^,_r."- 

7 l" — — 00 

Using the fact that I-n{z) = In{z), K-n{z) = Kn{z), we find that = T^. Hence, (1^ can be rewritten as 

OO 

M , —T” TT^'^ rpn"Tj21 

Mn,n' 2_^ '-^n,n"^2 n" ,n' 

n" — — OO 

oo 

=Tr ^ K^+^.TfK^,+r.". 


(32) 


(33) 


n — — OO 


III. LARGE SEPARATION ASYMPTOTIC BEHAVIOR 


In this section, we compute the asymptotic behavior of the Casimir interaction energy when L ^ R and L ^ i?i+i? 2 , 
for the cylinder-plate interaction and the cylinder-cylinder interaction when the two cylinders are outside each other. 
First notice that by making a change of variables 7 = 7 /L and expanding the logarithm, we have 

Ece. = ^ ^ In (1 - M) 

^0 


LD-i JO 

hcH 


^ s-k 1 


fi77 


D-2 


E E 

nQ — — oo Us—— 00 


(34) 


For the cylinder-plate case, 

= T^ftK^+n'{2^). 

For the case of two cylinders exterior to each other, 

OO OO 

^ A„+„.( 7 )rfx„,+„»( 7 )= ^ 


no,ni ■ • ■ ■ 


n" — — 00 


n"=—oo 


Here 


and 


rjinX> _ d-n i^Ri/E) 


Kn {lR^/L) 


^n,N _ i^Eti/L) 




= 7;"A„+„.(7). 

From these, we see that to determine the large separation asymptotic behavior (i.e. L ^ R in the cylinder-plate 
case and L ^ Ri + R 2 in the cylinder-cylinder case) of the Casimir interaction energy, we need to find the small z 
asymptotic behavior of /„(z)/A„(z) and I'^{z)/K'^{z). 

From any standard textbook of special functions, we find that 

4(^) 1 I 

Inz ’ 

1 o 


Ko{z) 


K'oiz) 2 ^ 


'■OV 

I[{z) 


1 


K[{z) 2 ^ 

In{z) 


Kn{z) 

KM. 

K'niz) 


=0 (z^") , n > 1, 

=0 (z^") , n > 1 . 


(35) 
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Hence, for the cylinder-plate case, we find that when the cylinder is imposed with Dirichlet boundary conditions, 
the leading term of the large separation asymptotic expansion comes from the term s = 0 and riQ = 0; whereas when 
the cylinder is imposed with the Neumann boundary conditions, the leading term of the large separation asymptotic 
expansion comes from the term s = 0 and ng = 0, ±1, namely. 




iDD 

Cas 


hcH 


pDN 

^Cas 


pND 

■^Cas 


pOO 


hcH 


E, 


-NN 

Cas 


2 D-in^r {E^) LD-i Jo V O'O ’ 

2D-i^^r(d^)LD-i Jo V O'O bi -I’-V 


Straightforward computation then gives 


E, 


DD 

Cas 


hcHT {E^) 


2 D+in—LD-nn{L/R) 


E, 


DN 

Cas 


hcHT{E^) 


pND 

^Cas 


2 D+i-jr—LD-Hn{L/R) 
hcH{3D + i)r (^) i?2 


2^+37r'^L^+^ 


E, 


-NN 

Cas 


hcH{^D + i)r (^) R 

2^+37r^^2^ L^+i 


(36) 


Notice that if Dirichlet boundary condition is imposed on the cylinder, the leading term is of order L“^+^/ln(L); 
whereas if Neumann boundary condition is imposed on the cylinder, the leading term is of order L~^~^. 

When D — 3, (I5S1) reads as 


pDD 

^Cas 


hcH 


T7' 

tj. 


NN 

Cas 


In (L/R) ’ 

3hcHR^ 

32'kL^ 


(37) 


These agree with the results obtained in 

For two cylinders that are exterior to each other, we rewrite (IMl) as 


Ecas — 


hcH 


2 ^-%^r(^)L^-i 


E 

s=0 


5 + 1 




D-2 


V ... V V ... V ...N^ . 

/ > / > / > / > ^’no,nQ ’nQ,ni ■ • ■ ^’ns,n+’n^,no 

no = — oo ns——cx>nL = — oo n'— — c<D 


(38) 


As in the cylinder-plate case, we find that the leading terms of the large separation asymptotic expansions are given 

by 


pDD 

^Cas 


piDN 

^Cas 


hcEl 


2 D-i^^Y LD-i 


7 D — 2 Arl.D Ar2,]D 

dll N),o N),o ’ 


hcEl 


f ^ — -D — 2 ( atIiD a7'2,N I 7\t 1,D Ar2,N . 7 \t1,D at2,N \ 

- d77 [Nofi ^ 0,0 + + d^o:-id^-\,o) , 


77 

tj, 


NN 

Cas 


2 ^-%—r(^)L^ 

/'°° ,__D_2 /,7,N,r2,N , A7,N,,r2,N , ,rl,N ,,r2,N 

■ 2^-l^^r(j^)L^-l Jo ^1.0 + 

+N^_foNo:-i + E.’o^^o.f + ++,f+ +A^+i+;!'i) , 

















and is obtained from by interchanging Ri and i? 2 - 
Straightforward computation gives 
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iDD 

hcHT{i^) 

(Jas 

2D+Itt^T (f) In (L/Ri) In {L/R 2 ) 

^DN 

hcT {3D + l)HRl 

'Cas 

2D+37r^r (^) L^+i In {L/Ri) ’ 

rNN 

3tic{D + 1){3D + 5)r {R^)'^HRlRl 

(Jas 

(^) LD+3 


(39) 


Notice that the leading term of the Casimir interaction is of order L“^+^/(lnL)^, L“^“^/(lnL) and L~^~^ respec¬ 
tively for DD, DN and NN boundary conditions. 

From above, we see that at large separation, the decay of the Casimir interaction is slower when Dirichlet boundary 
condition is imposed on the cylinder, and the decay is faster when Neumann boundary conditions is imposed on the 
cylinder. For the same boundary conditions, the decay is faster in higher dimensions. 


IV. SMALL SEPARATION ASYMPTOTIC BEHAVIOR 

Small separation asymptotic behavior of the Casimir interaction is of much more interest since Casimir force is 
inversely proportional to some power of the distance between the objects. It is always expected that the leading 
term of the Casimir interaction should agree with that derived using the proximity force approximation. A subject 
of much more interest is the next-to-leading order term, because the ratio of the next-to-leading order term to the 
leading order term is experimentally measurable. For the cylinder-plate configuration in (3 + 1) dimensions, the 
small separation asymptotic behavior has been derived in up to the next-to-leading order term. The idea of the 
derivation is to similar to perturbation in quantum field theory. Later the method has been generalized to compute 
the small separation asymptotic behavior of the Casimir interaction energy in various settings [^. 1^. 1^. 1545 Ij . 

In Section m we have seen that for the cylinder-plate or cylinder-cylinder configurations, the Casimir interaction 
energy is given by 


Acas = 


hcH 




d77-^"^Trln(I-: 


(40) 


with different matrix M for different scenarios. The first step in deriving the small separation asymptotic expansion 
is to expand the logarithm in (l40l) . which gives 


^ hcH ^ 1 

~ ~rD-i\ 71 


r ( -P-i j ^ s + l 
^ V 2 / s =0 


d'jj 


.-2 ^ 

no=—oo 


■ E 

ns — — co 


Afno.ni ■ • ■ ^Us ,nQ- 


(41) 


In the following, we will discuss the different scenarios separately. 


A. The cylinder-plate case 


In this case, define 


e 


d 

— , n = no. 


uj = yi?. 


and make a change of variables 


n-i = n + fii, 1 < i < s, 
nVl — 


UJ = 


T 
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Approximating the summation by corresponding integrations, we have 


^ hcH ^ 1 




dfll ... I dfls-^nQjTii ■ • • -^n 


where 


(2a;(l + £)). 

Here X = D or N is the boundary condition on the sphere, and Y = D or N is the boundary condition on the plate, 

Q !]3 = 0 , CTn = 1 , 


rpni,D _ Iriiiu}) 


rjirii.N _ 


Now we need to find the asymptotic expansion of Mrn^m+i keeping in mind that n~e hi~£ 2 . We also need 
the Debye asymptotic expansions of the modihed Bessel functions given below: 

Ii,{yz) ~ I 1 H-f • • •) 5 

\/2ttv (1 + V V ) 


where 


From this, we find that 


Therefore, 


where 


K^(yz) 


A uijtjz)) 




V 2i2(1 + z2)5 V / 

1 e'"''(^)(l + A)3 A t>i(f(z)) 

^/2ttv z \ V 


<(„). _ /Z£y2i!<i±£!)i f 1 _ nWfll +,. i 

\l 2 v z \ V ) 


ri{z) = a/I + + log 


1 +VT+A 


t{z) = 


, , t 5A 
“■''•= 8 “ 24- 


8 N 3f 7A 

”'<** = "s'+ 


AAA A ^ 2 m(t(z)) 

Ki,{yz) TT \ V 

AAA _ _ }_ 2^v{-) A , 2 ^^l(A^)) 

K[j{vz) TV \ V 


1 2!8i?7(zi)-l/2l7(22) 

MS.*. ~ (- 1 )”+” ^ (1 + Af). 

y'ZTTl22 (1 + Z2)4 


I 2 i = Ui, 122 = ni+ Ui+I 

(jj 2a;(l + e) 

2^1 = —, 2^2 =-^-, 

rii Tii + rii+i 


^2 =- (^ 2 ui(A - 
A2 =- (2vi{t) - hi (A 
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are of order s. With the help of a computer symbolic math package, we find that 




'Tii —rii^i _ 


2 y nn 


— (1 + Bi^i + Si, 2 ) exp- 


2en T{ni-fii+iY 


An 


(l + Af), 


where Si,i and Si ,2 are respectively terms of order y/e and £. Substitute into (l42ll . we find that 


pXY 

^Cas 


hcH 




-E 




xexp[-^^^i±^-^ 

V ^ i =0 


S + 1 
2 ' 


2®7r2 


dn 




,-D-l 


dr- 


D-i 


dfii ... 


dfis 


T{ni - fij+i) 

4n 


T 

s — 1 s 


(1 + (s + l)v4^) j 1 + Si,i + Si,iSj,i + Si ,2 

2^0 


2=0 


2=0 


The integration over hi is Gaussian, and it has been explained in (see also [i^). One finds that the terms of 
order y/e would not contribute since it is odd in one of the hi. After the integration, one is left with an expression of 
the form 


pXY 

^Cas 


hcH 




(_2((“x+aY)(s+l) 

(7TT)i 


r. The ini 
hcT {D - i) Hy/R 


dr 


T^-i 


dn 1 


exp 


2£(s + l)n 


(1 + J-^), 


where is a term of order e. The integration over n is straightforward using the definition of gamma function. One 
obtain 


pXY 


^ (« + iP+ 


nE 


(^_]^^(Q^x+aY)(s+l) /*! 


s+J-j r) . 

— dr{l-r^)^{l + g^) 

Jo 


where is a term of order e and is a polynomial of degree two in r^. Using 


r(efs) 


we find that 


where 


pXY 

^Cas 


her {D - 1) HVR - (-l)(-^x+.Y)(a+l) ^ 

22z.-i^^r (f) d^-i (s +1)^+^ ^ ^ 


s=0 


= 


^ ( 4^-5 _ (^- 2 )P- 3 ) 2 \ 

Vl2(2S-3) 3S(2S-3) ^ ’J 

/ 4S-5 D^ + 7D-6 d 

V 12(2S - 3) ~ 3D{2D JR 


d_ 

S’ 


Finally using the fact that 


E 

s=0 


(s + iy 


= CW. E77:rTV = -(i-2 


-fe+l'l 


s =0 


(s + 1)'' 


C(fc), 


we have 


pDD 

^Cas 


hcT {D - i) C(S + 1)hVR 




1 + 


4S-5 


(^-2)(S-3) C(i^-l) 


£;DN _ o her {D-^)aD + l)Hy/R / 

'' ’ ') 2 D-A.^^y'(R) d^~^ V 


12(2S-3) 3D{2D-3) C(-D + 1)JS 

4S-5 (S-2)(S-3)2^-4C(S-1) 


+ 


2 .ND n ^_^,her{D-^)aD + l)HVR 

^Cas ^ ) „,r> 1 D-i„ /n\ .n i ^ 


12(2S-3) 3D{2D-3) 2^-lC(S + l)JS "j" 

AD-5 + 7D-62^-AC{D-iy 


22 D-l^^j^(D)dD-^ 


(:D) V [l 2 ( 2 S - 3) 3D{2D - 3) 2 ^ - 1 C(S + 1)J R 


+ ... 


TA 

-C/, 


'NN 

Cas 


her {D - i) C(S + l)HyjR 


220-l7r^r(f)d^-i 


I— C{D + 1)1 + 


S2 + 7S-6C(S-1) 


4S- 5 

12 ( 2 S- 3) 3S(2S-3) C(S + 1)J S ' "V ' 


(46) 
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It is easy to check that the respective leading terms coincide with the result of proximity force approximation (see 
Section IVJ. Hence, we can write 


where 


E 


XY 

Cas 


E 


PFA,XY 

Cas 





pPFA.DD _ pPFA.NN 
-^Cas — ^Cas 


pPFA.DN _ pPFA,ND 
^Cas — ^Cas 


her (D - i) c(E + 

22 D-i^£i^r (f) ’ 

_ ncr (D - ^) c(D + i)hVr 


^XY ^ ^Y^ 


^DD _ 4D — 5 
12(2D- 3) 
_ 4L> — 5 
12(2D- 3) 
^nd _ 4Z1 — 5 
^ “l2(2D-3) 

NN _ 411 — 5 

12(2D- 3) 


{D-2){D-3) gP-l) 
3D{2D-3) CP + 1)’ 

{D - 2){D - 3) 2^ - 4 C{D - 1) 
3D{2D-3) 2^-lC(£> + l)’ 
D‘^ + 7D-6 2° -4:C{D- 1) 
3D{2D-3) 2° - 1C(£> + 1)’ 
D^ + ID-Q C{D - 1 ) 
3D(2D-3) CiD + l)' 


(47) 


The values of are tabulated in Table U of Appendix El for 3 < D < 6. d measures the correction to the proximity 
force approximation. The dependence of i? on H is plotted in Fig. [T] 
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FIG. 1: The dependence of d on dimension D for different combinations of boundary conditions. 


One observes some interesting phenomena, d is positive when the cylinder is imposed with Dirichlet boundary 
conditions, which indicates a positive correction to the proximity force approximation; and d is negative when the 
cylinder is imposed with Neumann boundary conditions, which indicates a negative correction. The correction is 
larger in the latter case. However, for all combinations of boundary conditions, we find that the magnitude of the 
correction decreases with dimension D. In fact, from (H71) . we find that when H ^ 1, 

^NY__111 (48) 

8 D’ 8 D’ ^ ^ 

which is inversely proportional to D. This is in big contrast to the sphere-plate interaction [4l|, where it is found that 
d ^ —D/A regardless of the boundary conditions, and thus becomes negative when D is large. In the sphere-plate 
case, the correction to the proximity force approximation becomes large when the dimension of spacetime is increased. 
However, for the cylinder-plate case, the correction to the proximity force approximation becomes smaller when the 
dimension of spacetime is increased. 
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B. The case where one cylinder lies parallelly inside the other 


In this case, 




Define 


and make a change of variables 


n-— — oo 


R 




R 2 


e = 


i ?2 — Ri R2 ~ Ri R2 ~ Ri 

n = no, LO = 7(i?2 - -Ri), 


m = n + fii, 1 < i < s, 

, b b ~ ^ 

ni = -n+ — (jii + Ui+i ) + Qi, 
a 2a 

nVl — 


or 


Approximating summations by integrations, we find that 

Ea„ - 

2 D- 27 r—S+1 Jo Jo 


D 


dTl\ ... I A /720 ,ni ■ • ■ -^n 


where 


, E (^(1 - ^)) . (^(1 - ^)) > 


rL.— — oo 


rj^ni,D Ri 

^ Km (aoj) ’ 

fnlB Kn'.ibuj) 

2 


' K'miau^y 

_ KJbu.) 

Ir, - 


lyibuj) ■ 

Using Debye asymptotic behavior of modified Bessel functions (H51) . we find that 

2 g2!/ir)(zi) + l/2»?(z2)-2!^3';(z3) + '^4»)(z4) 


r° 

J — c 


dqi 


2 ' Ky / v 2 i’i (1 + 21)4(1 + 24)4 


(1 + A2 + ISy + C2) , 


where 


vi=ni, i 22 =ni-n„ V 3 = n„ V4 = ni-n,+i, 

out a;(l — e) btu a;(l —e) 

zi = —, 2:2 = -, 23 = —, 24 = -; 

I'l V2 ^'3 Vi 


A 2 =-Mi(t), 
n 

C 2 = - Ui(T). 

n 


•^2 = 

n 


■s 5^0 5 

(49) 


(50) 


(51) 


(52) 
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A 2 , B 2 and C 2 are terms of order e. As in the cylinder-plate case, expanding each term keeping in mind that n has 
order and fii and qi has order , we obtain 


/ OO / 

dft^(l + iPz.i+iP*.2)expf- 
X {1 + Af +Bl +C 2 ), 


2 en br , 

-- 1 — 

ar 4n bn 


(53) 


where T>i i and T>i ^2 are respectively terms of order -y/e and e. The integration over qi is straightforward and gives an 
expansion of the form 




»2(_i)ax+«Y (1 -I- Oi l + t/i 2)exp {l + A 2 + B^ +C 2 ) . (54) 

2 v 7 rn \ or 4n / 


The rest is similar to the cylinder-plate case. We find that the up to the next-to-leading order term, the Casimir 
interaction energy can be written as 


77’DD _pDD,PFA 

-^Cas —-^Cas 


pDN _TpDN,PFA 

^Cas —^Cas 


pND _pND.PFA 

^Cas —^Cas 


pNN _pNN,PFA 
^Cas -^Cas 


1 -b 


1 + 


1 + 


1 + 


4£»- 5 


4(2L> - 3) i?2 - i?i 
AD-b d 


;^ddA_^dd d 


Ri 


R 2 

d 


4(2£> - 3) A 2 - Ri Ri i?2 y ’ 

AD — 5 d . wn d dn d 


A{2D - 3) A 2 - Ai 




4£»- 5 


i?2 

d 


- - - + — 

A{ 2 D - 3) i?2 - Ri Ri R 2 


Here, is defined in (HTl) . and are equal to the for the cylinder-plate case, and ^ is the leading term 

that coincides with the proximity force approximation. They are given explicitly by 

hcHT {D - ^) CiD + 1) I R 1 R 2 


E, 


_ 77’ 

= -C/, 


NN 


DD _ 

Cas,PFA — -^Cas,PFA 


pDN _ Z^ND _/-I 9 

'^Cas,PFA — '^Cas,PFA —U “ ^ 


220-W^r (f) d ^-5 V i ?2 - i?l ’ 

_ hcHT{D-^)aD + l) I 


2 ^^ 2 7r^2 (-y) d^ 2 \ R 2 — Rl 


Hence, we have 


pXY 

^Cas 


E, 


PFA.XY 

Cas 


= n -H d 


XY 


i ?2 — Rl V R 2 ~ Rl 


where 




4£»-5 
A{2D - 3) 


^ddI _^dd1 
h 


^DN ^ 414 5^ + ;<DNi _ 1 

b 




4(214 - 3) 
414-5 
4(214 - 3) 


a 

1 

a 


,NDi_^Nl 


(55) 


nn _ 414 5 1 _ 1 

4(214-3) a b' 

Recall that 5 = a -I- 1. Hence, we can regard d as depending on dimension D and /3 = 6 /a = l? 2 /l?i-the ratio of 
the radius of the larger cylinder to the radius of the smaller cylinder. Then 


1 


- 1 


6 = 


'- 1 


In Fig. [21 we plot the dependence of d on the dimension D and radii ratio jd for different boundary conditions. 
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FIG. 2: The dependence of i9 on dimension D and radii ratio /3 for different combinations of boundary conditions. 

We observe that i? is always positive when the inner cylinder is imposed with Dirichlet boundary conditions. When 
the inner cylinder is imposed with Neumann boundary conditions, t? can be positive or negative depending on the 
dimension D and the ratio of the radii of the cylinders. When D is large, we observe some universal behavior. In 
fact, from (l55ll and the asymptotic behavior of obtained in (j48]), we find that when 79 ^ 1, 




XY 


1 

' 2 ' 


regardless of the boundary conditions. This agrees with the graphs we obtained in Figure [2] The dominating term 
actually comes from 


479-5 
4(279 - 3)’ 


which is universal for all boundary conditions. 


C. The case where two parallel cylinders are exterior to each other 


In this case, 




n — — GO 


Define 


R 




R2 


e = 


R\ H“ R2 Ri R2 Ri R2 

n = no, io = 7(7 ?i + R 2 ), 
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and make a change of variables 


TLi = n + fii, 1 < i < s, 

b b 

-n+ — 
a 2 a 


, b 5 - N 

n-i = -n+ — (jii + m+i) + qi, 


n 


VT^- 


As in the previous case, we find that 
A^Cas 


HcH ^ 1 

s + 1 




(jJ = 


dn 'nP ^ dr 




dfi\ . . . / yjp, 

) P —CO 

(56) 


where 


= Tp ’^ E (^(1 + ^)) (^(1 + £)), 


n. — — co 


^ Kui {auj) ’ 

/yiTli ,N 

±1 — 

,^n',D 4'(&w) 

2 -Rr^ybu^y 

-^2 — 


K'^yau^y 

I'n'ibu;) 


KyXbu^)- 

Using Debye asymptotic behavior of modified Bessel functions (l4^ . we find that 


^(_l)«x+av / 

J —( 


dqi 


2 ^ ^2vir){z-L)-V2’ri{z2)+2v3ri{z3)-V4,ri{z4,) 

2tt 3/122124, (1 + z|)i(l + z|)i 


(1 + A 2 + By + C 2 ) , 


where 


i2i = Hi, j/2 = n'+nj, 123 =ni, i24 = p + n^+i, 


zi = 


au> 


vi 


Z2 = 


V 2 


Z3 = 


buj 

V 3 


Z4 = 


E±£). 

124 


(57) 


(58) 


A/ =-ui(r), 
n 

^?=^.r(r), 

r f \ 

€2 = ui{t). 

n 


n 

B, = ^vAt). 


As before, expanding each term according to orders of e gives 


Mr,. 


fjn ^2(•_2)ax+aY / (1 _)_p. 2)exp (-—--^(rii - Ui+i)^- 


' 27m 


2 en br 


ar An 


a T 
bn 


X {l + Ay +By +C 2 ) 


(59) 


(60) 


where Vi^i and are respectively terms of order r/e and e. The rest is similar to the case where one cylinder is 
inside the other. We find that the up to next-to-leading order term, the Casimir interaction energy can be written as 


pDO _pDD,PFA 

-^Cas -^Cas 


pON _pDN,PFA 

^Cas -^Cas 


pND _pND,PFA 

^Cas —^Cas 


pNN _pNN,PFA 
^Cas —^Cas 


1 - 


1 - 


1 - 


1 - 


5 d J3J3 /y , DD _y_ 

4(2B - 3) i?i + i?2 i?i R 2 J ’ 

—-— + 

4(2D - 3) i?i + i?2 i?i E 2 J 

— 5 _ d . ,dn ^ \ 

4(2D - 3) i?i + i?2 i?i R 2 )' 

AD — 5 d ,jvjN ^ I _y_\ 

A{2D - 3) i?i + i?2 i?i R 2 )' 


(61) 
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XY PFA 

Here, is the leading term that coincides with the proximity force approximation. They are given explicitly 


by 


Hence, we have 


where 


pDD _ pNN 

^Cas,PFA — -^Cas.PFA 


hcHT {D - b) ({D + 1) / R 1 R 2 


E, 


DN 

Cas.PFA 


= E, 


ND 

Cas.PFA 


=( 1-2 


22D-i^^^r(f )y R 1 + R 2 ' 

hcHT {D - i) CiD + 1) / R 1 R 2 


r2D-i 


\ R 1 +R 2 


pXY 

^Cas 


E, 


iPFA.XY 

Cas 


= n +1? 


XY 




^DN^_ 


Rl + i?2 

4L>-5 


Rl + R 2 


A{2D - 3) 
4L>-5 


+ ^Di + .,DD 


a 0 

+ ^n1 ^d1 

4(244-3) ^ a 6’ 


(62) 
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FIG. 3: The dependence of 2? on dimension D and radii ratio A for different combinations of boundary conditions. 


Recall that & = 1 — a. Hence, we can regard d as depending on dimension D and A = b/a = i?2/Ri“the ratio of 
the radii of the cylinders. Without loss of generality, we can assume that R 2 > Ri- Then A > 1, 

A 


1 


A + r 


b = 


A + l 


a = 
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In Fig. |3l we plot the dependence of i? on the dimension D and radii ratio A for different boundary conditions. 

We observe that d is always negative for DN, ND and NN boundary conditions. When both cylinders are imposed 
with Dirichlet boundary conditions, d can be positive or negative depending on the dimension D and the ratio of the 
radii of the cylinders. When D is large, we observe some universal behavior. In fact, from (1621) and the asymptotic 
behavior of obtained in (l48l) . we find that when ^ 1, 


regardless of the boundary conditions. This agrees with the graphs we obtained in Figure [H The dominating term 
actually comes from 


AD-b 

A{2D-?,y 


which is universal for all boundary conditions. 


V. POSTULATE FOR DERIVATIVE EXPANSION FORMULA 

In a series of papers [12, [13, HE] , Fosco, Lombardo and Mazzitelli used derivative expansion to compute the Casimir 
interaction energy between a curved surface and a plane with Dirichlet boundary conditions up to the next-to-leading 
order term. In (52l| . they showed that the derivative expansion of the Casimir interaction energy is given by 

Ey.1 = hc[ ^ + 62(D) (63) 

Js \ |^(xu)| l^(x±)| / 

Here x_l = {x 2 , ■ ■ ■ ,xo), Xi = '(/'(xj_), xj^ S S defines the position of the curved surface with respect to the plane at 
Xi = 0, and 


60(D) = - 

h 2 {D)=- 


r(i^)C(D + i) 


3 X 2^+'^'K^y 


(^) an - 1) + (i. + i)r {^1 1 az, +1) 


=6o(D) 


D + 1 (D-2)(D-3)C(D-I) 


6 6D C(D +1) / ■ 

Consider the cylinder-plate interaction. We can take the cylinder to be (xi — L)^ + x\ = B?. Then 

■!/’(xj_) ='ip{x 2 ) = L - \JR'^ - xl, 


Vtp = 


X2 




?e2- 


Some computations give 

f d^-ixi- 


rR 


IV'(x_l)| 


D 


^2H / dx2- 


2H^/R ttF {D - i) 

= 2^-iF(f)r(AM) V" ^ AR2D-3 " '' 7 ’ 


1 - 


3d 1 


(64) 


jD-1 


1 


rR 




IV’(xj.)| 


D 


=2H 


dX 2 


(l - 7i?2 - xi) 


D 


Rl-xl 


2Hy/R ttF (D - i) 2 d ^ 

72d^-5 2^-iF(f)r(^) 2D-3D^---' 
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Hence, derivative expansion gives 


E, 


DE,DD 

Cas 


hcT (D - i) ({D + 1)hVR { 


AD-b 


{D-2){D-3) CiD-l) 


12{2D-3) 3D{2D-3) ({D + 1) _ 


d 

R 


(65) 


agreeing with the result we obtain in the first formula of (H51) . 

Encouraged by this, we would like to give a postulate for the result of derivative expansion for DN, ND and NN 
boundary conditions, in the case where the two interacting objects are both curved. Inspired by [H^, let us formulate 
the following ansatz for the small separation asymptotic behavior of the Casimir interaction energy between two 
curved objects in {D + l)-dimensional Minkowski spacetime: 


(l + Pi(H)VHi ■ Vili + P2{H)VH2 ■ VH2 + /3x ■VH 2 + ..), (66) 


where is the Casimir energy density between two parallel plates, E can be taken to be the xi = 0 plane 
parametrized by xj^ = {x 2 , ■ ■ ■ ,xd), xi = Hi(xj^) and X 2 = H 2 (xj_) are the height profiles of the two objects with 
respect to E, and H = Hi — H 2 is the height difference. When the second object is a plane we can take it as the 
plane a;i = 0 and then H 2 = 0. 

Notice that 




"'Cas 


{H) = 


HD 


(67) 


where 


bDD^D) = briD)=- 


HcT {d^) CjD + 1) 




bnD)^briD)={l-2-^) 


Her (D^) ({D + 1) 


( 68 ) 


2£>+i7r^ 


According to the result of 


we mentioned above, 

D + 1 {D-2){D-3)aD-l) 


oDD 

^1 


6 


6 D 


CiD + l)- 


(69) 


We are going to determine the values of Pi, P 2 and /3x for DD, DN, ND and NN boundary conditions based on 
our results on the small separation asymptotic expansions of the Casimir interaction energies between two cylinders 
exterior to each other. 

We can take the two cylinders to have height profiles 


Hi =Li- jRf-x: 


H 2 — — L 2 + \ R 2 — x^. 


Then L = Li + L 2 \s the distance between the centers of the cylinders. Assuming that Ri < R 2 , some tedious 
computations give 


^D-l, 


x± 


1 

llD 


=2H 


H 


dx- 


1 


(^L — ijR\ — x"^ — \/R2 — x'^'^ 


D 


2 -kRiR 2 r (D - i) 
V i?i + i ?2 r(D) 


(70) 


1 + 


4(27? - 3) i?i + i?2 4(2D - 3) \Ri R 2 


,D-i, V77i • VHi 


rRi 


X±- 


HD 


=2H 


H 


dx- 


(l - ^Rj - x^ - ijRl - 
27ri?ii?2 r (D - i) 2 


D r2 


Ri-x^ 


dD-hy R1+R2 r{D) 2 D- 3 \Ri R1+R2 


+ .. 


(71) 
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rRi 


RD 


dx- 


RI-: 


H 




(l - ^Rj - - yjRl - x2) 

27ri?ii?2 r (£> - I) 2 / d d 

d?l + i?2 r(Zl) 2Z1 — 3 \i?2 + ^2 


+ .. 


(72) 




x±- 


Vidi • Vid2 


id^ 


= -2H 


dx- 


(l _ _ ^2 _ ^^2 _ 2 . 2 ^ ^Rl-X^^Rl-X^ 

H /27ri?ii?2 r(ll- i) 2 


(73) 


d^-5Vi?i+i?2 r(£l) 2L»-3Vi?i+^2 

Substituting into (I66L we find that for the interaction of two cylinders, 


pDE,XY _,XY 
^Cas —"0 


id 


27rd7id72 r (d? - i) 


d^-Hd?i + d ?2 r(dd) 

flXY/ 


1 + 


9 


4(2d? - 3) d?i + d?2 4(2D-3)V^i -^2 


+ ... 


d d 
+ 


(d 

d ^ 

1 , 2 yf^{D) , 

(d 

d \ 

2 PViD) d \ 

\Ri 

Ri + R 2 ) 

1+ 2D-3 ' 

\Ri 

Ri + R 2 J 

2D-3 d7i-fd?2 ■■■] 


= 6 , 


XY 


id 


27rd7id72 r(d?-i) 1 


+ - 


d^-5Vi?i + i ?2 r(D) 

1 


2dd - 3 V 4 ■ {D) - 2 PYiD) 


(74) 


i?i H" -H 2 




2D-3 

Compare to our results (EB), we find that the leading terms do agree, and the next-to-leading order terms give 
j - 213^(0) - mriD) - 213"(D) = 

Y1 ^ SD + \ ) Y1 V ;> 

DN 3_4L»-5 (id-2)(d?-3)C(7?-l)2^-4 


2/3)^^^ - - = 

- I = 

o oNN 3 _ 
^Pl 1 - 


12 312 + 22.31' />r(Z2) = ,2r(«, 

4d3-5 £)2-H7dl-6C(7?-l)2^-4 


(75) 


12 


3D 


C(dl-f 1)2^-1 


, /3rp) = /3r(id) 


4dl - 5 D"^ + 1D -&C{D - 1) 

3d3 C(i^ + 1)' 


/32^^(dd)=/3r ( 2 D). 


From these, we obtain 


and 


/^^(^D) = - ^^(.D) - /3^^(D), 


(76) 


/3?°(id) =/32°°(dd) = 
/3?N(dd) = /32ND(^) ^ 

/3r(2D) = ^F^(D) = 

;9r(2D) = ;92^''(dd) = 


£»-fl (D - 2)(D - 3) ((D - 1) 

6 6 D ((D + 1)’ 

D + 1 (D-2)(D-3)((D-l)2^-4 

~6 6id C(dl-f 1)2^-1’ 

D + 1 D^ + 7D - 6 ((D - 1) 2^ - 4 
~6 67d C(2D-f 1) 2^ - 1 ’ 

D + 1 D^ + 7D - 6 ((D - 1) 


(77) 


6 6dl 

When D = 3, these agree with the results obtained in [s^. 


({D + iy 
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Now let us compare these results to the results of two spheres we obtained in [d^. In this case, we take the height 
profiles of the spheres with radii i?i and i ?2 to be 


i?i = Li - \ Rl - 




H 2 = -L 2 + \ Rl- x]_, 


where xi_ = + + xj^. Li + L 2 = L is the distance between the centers of the spheres. 

Now, assuming i?i < R 2 , some tedious computations give 


D—1 D 

1 27r— 




D-2 


(l - - xl - ^Ri - xi) 


D 


' 2 ^r{R)d^ \Ri + R 2 


f R 1 R 2 3(D + 1) d D + 1 / d d 

4 i?i + i?2 4 “ 


(78) 


jD-1 ViJi • ViJi 27r 2 


rRi 


d^-^K±- 


hd 


dx±x 


D-2 




D 

7T 2 


R 1 R 2 


2 ^r{R)d^ \Ri + R 2 j \Ri Ri + R‘. 


/ 2 d 2 d 


(79) 


,D-1 Vi72 • Vila _ Stt^ 




D-2 


(^L - ^/Rl-x\ - 


L2 u.j_ 


TT 2 




2 ^r(§)d^ \Ri + R2j \R2 R1 + R2 


2 d 2 d 


(80) 


ji’-i 


XX- 


Viii • VH 2 

ID 


27r-V 

W) 


rRi 


dxj_x^ 


D-2 


(h-^Ri-xi-^Ri-xiY ^Hi-xi^Rl-x^ 


TT 2 


R 1 R 2 


2^r(f)d^ \Ri + R2 

Hence, for the interaction of two spheres, our ansatz (l66l) gives 

3(11 + 1) d 


2 d 


Rl + i?2 


+ ... . 


c,DE,XY _;.XY 
-^Cas —"0 . 0-1 


2 ^T{R)d^ 

d 


+Dr{D) 


1 + 

d 


— h 

— L/i 


TT 2 


XY_ 

0 


2 ^^V (—) d^ 


Rl Rl + i?2 

1 + 


ii + 1 
4 i?i + i?2 4 

XY / f d d 


1 


+ 2/?2^Mii) 


i?2 iil + i?2 


-DriD) 


_ 2/3fY(7^) _ _ 2l3YiD) 


Rl + i?2 
d 


Rl + i?2 


+ 2/3f^(ii)- 


ii + 1 


— +[2priD)- 


D + 1 


d 


(81) 


(82) 


With the values of /3i and P 2 given by dzzi), this agrees perfectly with the result we obtained in for two spheres 
when D ^ A. 

Now some explanations are in order. The derivative expansion technique is a formal and non-rigorous method to 
obtain the small separation asymptotic behavior of the Casimir interaction energy. The result might not be correct 
due to some un-observed singularities in the formal derivation. Therefore, the ansatz (I66p can only be used as a 
reference for the small separation asymptotic expansion of the Casimir interaction energy, but it needs to be checked 
against actual computations. 
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VI. CONCLUSION 

In this work, we have considered the Casimir interaction in [D + l)-diinensional spacetime due to the vacuum 
fluctuations of massless scalar fields between a cylinder and a plate, between two parallel cylinders where one is inside 
the other, and between two parallel cylinders exterior to each other. We derive the explicit integral representations 
for the Casimir interaction energies and use them to study the large separation and small separation asymptotic 
behaviors of the Casimir interactions. The large separation asymptotic behaviors are easy to compute and the order 
of decay is smallest in the Dirichlet-Dirichlet case, and largest in the Neumann-Neumann case. The computations 
of the small separation asymptotic behaviors are more complicated. The leading terms are found to agree with the 
proximity force approximation. The results on the next-to-leading order terms are important and they exhibit some 
universal behaviors. In particular, we find that for the cylinder-plate case, the ratio of the next-to-leading order term 
to the leading order term is inversely proportional to D. For the case where one cylinder is inside the other, the ratio 
of the next-to-leading order term to the leading order term approaches the limiting value 1/2 when D is large. For the 
case where the two cylinders are outside each other, the ratio of the next-to-leading order term to the leading order 
term approaches the limiting value —1/2 when D is large. Hence, we find that the ratio is bounded in dimensions for 
all cases we consider. Therefore, the corrections to the proximity force approximation will not gets larger in higher 
dimensions, in contrast to the sphere-plate and sphere-sphere interactions, where it is found that the ratio of the 
next-to-leading order term to the leading order term is proportional to D when D is large [4l|, . 

An interesting thing to note is that our small separation asymptotic expansion for the case of Dirichlet-Dirichlet 
cylinder-plate interaction agrees with the result derived using derivative expansion in . Generalizing the D — 3 case 
in [s^ , we postulate a general form of the derivative expansion for small separation asymptotic expansion of the scalar 
Casimir interaction energy in (D + l)-dimensional Minkowski spacetime, for two curved surfaces with combinations 
of Dirichlet and Neumann boundary conditions, based on our results on the cylinder-cylinder interaction. We also 
check our postulate with the results we obtained for the sphere-sphere interaction in and find that the postulate 
gives correct expansion except when D = 4. 
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Appendix A: Tabulation of constants 


TABLE I: The values of for 3 < D < 6. 


D 


^DN 

^ND 

^NN 

exact 

numerical 

exact 

numerical 

exact 

numerical 

exact 

numerical 


7 


7 


7 160 


7 40 


3 


0.1944 


0.1944 


-0.5775 


-1.1565 


36 


36 


36 2l7r2 


36 37r2 


4 

11 1 C(3) 

0.1447 

11 2 C(3) 

0.1524 

11 38 C(3) 

-0.4040 

11 19C(3) 

-0.5509 


60 30 C(5) 


60 75 C(5) 


60 75 C(5) 


60 30 C(5) 



5 3 


5 84 


5 756 


5 27 


c: 


0.1178 


0.1237 


-0.3156 


-0.3686 


28 57r2 

28 1557r2 

28 1557r2 

28 57r2 

6 

19 2 C(5) 

0.0998 

19 40 C(5) 

0.1034 

19 80 C(5) 

-0.2594 

19 4C(5) 

-0.2811 


108 27 C(7) 


108 567 C(7) 


108 189 C(7) 


108 9 C(7) 
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